Continued improvement in the development of the SARS-CoV-2 whole genome sequencing proficiency testing program.

Application of whole genome sequencing (WGS) has allowed monitoring of the emergence of variants of concern (VOC) of severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) globally. Genomic investigation of emerging variants and surveillance of clinical progress has reduced the public health impact of infection during the COVID-19 pandemic. These steps required developing and implementing a proficiency testing program (PTP), as WGS has been incorporated into routine reference laboratory practice. In this study, we describe how the PTP evaluated the capacity and capability of one New Zealand and 14 Australian public health laboratories to perform WGS of SARS-CoV-2 in 2022. The participants' performances in characterising a specimen panel of known SARS-CoV-2 isolates in the PTP were assessed based on: (1) genome coverage, (2) Pango lineage, and (3) sequence quality, with the choice of assessment metrics refined based on a previously reported assessment conducted in 2021. The participants' performances in 2021 and 2022 were also compared after reassessing the 2021 results using the more stringent metrics adopted in 2022. We found that more participants would have failed the 2021 assessment for all survey samples and a significantly higher fail rate per sample in 2021 compared to 2022. This study highlights the importance of choosing appropriate performance metrics to reflect better the laboratories' capacity to perform SARS-CoV-2 WGS, as was done in the 2022 PTP. It also displays the need for a PTP for WGS of SARS-CoV-2 to be available to public health laboratories ongoing, with continuous refinements in the design and provision of the PTP to account for the dynamic nature of the COVID-19 pandemic as SARS-CoV-2 continues to evolve.

Human cytomegalovirus (CMV) dysregulates neurodevelopmental pathways in cerebral organoids.

Human cytomegalovirus (CMV) infection is the leading non-genetic aetiology of congenital malformation in developed countries, causing significant fetal neurological injury. This study investigated potential CMV pathogenetic mechanisms of fetal neural malformation using in vitro human cerebral organoids. Cerebral organoids were permissive to CMV replication, and infection dysregulated cellular pluripotency and differentiation pathways. Aberrant expression of dual-specificity tyrosine phosphorylation-regulated kinases (DYRK), sonic hedgehog (SHH), pluripotency, neurodegeneration, axon guidance, hippo signalling and dopaminergic synapse pathways were observed in CMV-infected organoids using immunofluorescence and RNA-sequencing. Infection with CMV resulted in dysregulation of 236 Autism Spectrum Disorder (ASD)-related genes (p = 1.57E-05) and pathways. This notable observation suggests potential links between congenital CMV infection and ASD. Using DisGeNET databases, 103 diseases related to neural malformation or mental disorders were enriched in CMV-infected organoids. Cytomegalovirus infection-related dysregulation of key cerebral cellular pathways potentially provides important, modifiable pathogenetic mechanisms for congenital CMV-induced neural malformation and ASD.

Shrinking in the dark: Parallel endosymbiont genome erosions are associated with repeated host transitions to an underground life.

Microbial symbioses have had profound impacts on the evolution of animals. Conversely, changes in host biology may impact the evolutionary trajectory of symbionts themselves. Blattabacterium cuenoti is present in almost all cockroach species and enables hosts to subsist on a nutrient-poor diet. To investigate if host biology has impacted Blattabacterium at the genomic level, we sequenced and analyzed 25 genomes from Australian soil-burrowing cockroaches (Blaberidae: Panesthiinae), which have undergone at least seven separate subterranean, subsocial transitions from above-ground, wood-feeding ancestors. We find at least three independent instances of genome erosion have occurred in Blattabacterium strains exclusive to Australian soil-burrowing cockroaches. These shrinkages have involved the repeated inactivation of genes involved in amino acid biosynthesis and nitrogen recycling, the core role of Blattabacterium in the host-symbiont relationship. The most drastic of these erosions have occurred in hosts thought to have transitioned underground the earliest relative to other lineages, further suggestive of a link between gene loss in Blattabacterium and the burrowing behavior of hosts. As Blattabacterium is unable to fulfill its core function in certain host lineages, these findings suggest soil-burrowing cockroaches must acquire these nutrients from novel sources. Our study represents one of the first cases, to our knowledge, of parallel host adaptations leading to concomitant parallelism in their mutualistic symbionts, further underscoring the intimate relationship between these two partners.

Clinical, Genomic, and Immunological Characterization of RSV Surge in Sydney, Australia, 2022.

OBJECTIVES: The 2022 seasonal respiratory syncytial virus (RSV) epidemic in Sydney, Australia saw an unprecedented number of RSV detections. We aimed to characterize genomic and immunologic factors associated with the surge in RSV cases. METHODS: Whole genome sequences of RSV were generated from 264 RSV-infected infants and linked to case-matched clinical data from the 2022 southern hemisphere RSV season. We then performed an immunologic analysis of baseline RSV-specific humoral immunity in women of childbearing age before and throughout the coronavirus disease 2019 pandemic. RESULTS: Clinical analysis revealed a high burden of disease across patients of all health backgrounds. More than one-half of RSV-related health care visits by infants resulted in hospitalization, and one-quarter required high-flow respiratory support or a higher level of care. Viral phylogenetic analyses revealed that 2022 Sydney RSV sequences were closely related to viruses that had been circulating globally since 2017, including those detected in recent US outbreaks. Nonsynonymous mutations within the palivizumab and nirsevimab binding sites were detected at low frequencies. There was no difference in baseline RSV-neutralizing antibody titers between 2020 and 2022. CONCLUSIONS: Collectively, these findings suggest that neither the emergence of a novel RSV genotype nor hypothesized immune debt was associated with the surge of RSV cases and hospitalizations in 2022. Continued genomic and immunologic surveillance is required to further understand the factors driving outbreaks of RSV globally, and to inform guidelines for the rollout and ongoing use of recently developed immunotherapeutics and vaccines.

Bioinformatic investigation of discordant sequence data for SARS-CoV-2: insights for robust genomic analysis during pandemic surveillance.

The COVID-19 pandemic has necessitated the rapid development and implementation of whole-genome sequencing (WGS) and bioinformatic methods for managing the pandemic. However, variability in methods and capabilities between laboratories has posed challenges in ensuring data accuracy. A national working group comprising 18 laboratory scientists and bioinformaticians from Australia and New Zealand was formed to improve data concordance across public health laboratories (PHLs). One effort, presented in this study, sought to understand the impact of the methodology on consensus genome concordance and interpretation. SARS-CoV-2 WGS proficiency testing programme (PTP) data were retrospectively obtained from the 2021 Royal College of Pathologists of Australasia Quality Assurance Programmes (RCPAQAP), which included 11 participating Australian laboratories. The submitted consensus genomes and reads from eight contrived specimens were investigated, focusing on discordant sequence data and findings were presented to the working group to inform best practices. Despite using a variety of laboratory and bioinformatic methods for SARS-CoV-2 WGS, participants largely produced concordant genomes. Two participants returned five discordant sites in a high-Cτ replicate, which could be resolved with reasonable bioinformatic quality thresholds. We noted ten discrepancies in genome assessment that arose from nucleotide heterogeneity at three different sites in three cell-culture-derived control specimens. While these sites were ultimately accurate after considering the participants' bioinformatic parameters, it presented an interesting challenge for developing standards to account for intrahost single nucleotide variation (iSNV). Observed differences had little to no impact on key surveillance metrics, lineage assignment and phylogenetic clustering, while genome coverage <90 % affected both. We recommend PHLs bioinformatically generate two consensus genomes with and without ambiguity thresholds for quality control and downstream analysis, respectively, and adhere to a minimum 90 % genome coverage threshold for inclusion in surveillance interpretations. We also suggest additional PTP assessment criteria, including primer efficiency, detection of iSNVs and minimum genome coverage of 90 %. This study underscores the importance of multidisciplinary national working groups in informing guidelines in real time for bioinformatic quality acceptance criteria. It demonstrates the potential for enhancing public health responses through improved data concordance and quality control in SARS-CoV-2 genomic analysis during pandemic surveillance.

Emergence and antibody evasion of BQ, BA.2.75 and SARS-CoV-2 recombinant sub-lineages in the face of maturing antibody breadth at the population level.

BACKGROUND: The Omicron era of the COVID-19 pandemic commenced at the beginning of 2022 and whilst it started with primarily BA.1, it was latter dominated by BA.2 and the related sub-lineage BA.5. Following resolution of the global BA.5 wave, a diverse grouping of Omicron sub-lineages emerged derived from BA.2, BA.5 and recombinants thereof. Whilst emerging from distinct lineages, all shared similar changes in the Spike glycoprotein affording them an outgrowth advantage through evasion of neutralising antibodies. METHODS: Over the course of 2022, we monitored the potency and breadth of antibody neutralization responses to many emerging variants in the Australian community at three levels: (i) we tracked over 420,000 U.S. plasma donors over time through various vaccine booster roll outs and Omicron waves using sequentially collected IgG pools; (ii) we mapped the antibody response in individuals using blood from stringently curated vaccine and convalescent cohorts. (iii) finally we determine the in vitro efficacy of clinically approved therapies Evusheld and Sotrovimab. FINDINGS: In pooled IgG samples, we observed the maturation of neutralization breadth to Omicron variants over time through continuing vaccine and infection waves. Importantly, in many cases, we observed increased antibody breadth to variants that were yet to be in circulation. Determination of viral neutralization at the cohort level supported equivalent coverage across prior and emerging variants with isolates BQ.1.1, XBB.1, BR.2.1 and XBF the most evasive. Further, these emerging variants were resistant to Evusheld, whilst increasing neutralization resistance to Sotrovimab was restricted to BQ.1.1 and XBF. We conclude at this current point in time that dominant variants can evade antibodies at levels equivalent to their most evasive lineage counterparts but sustain an entry phenotype that continues to promote an additional outgrowth advantage. In Australia, BR.2.1 and XBF share this phenotype and, in contrast to global variants, are uniquely dominant in this region in the later months of 2022. INTERPRETATION: Whilst the appearance of a diverse range of omicron lineages has led to primary or partial resistance to clinically approved monoclonal antibodies, the maturation of the antibody response across both cohorts and a large donor pools importantly observes increasing breadth in the antibody neutralisation responses over time with a trajectory that covers both current and known emerging variants. FUNDING: This work was primarily supported by Australian Medical Foundation research grants MRF2005760 (SGT, GM & WDR), Medical Research Future Fund Antiviral Development Call grant (WDR), the New South Wales Health COVID-19 Research Grants Round 2 (SGT & FB) and the NSW Vaccine Infection and Immunology Collaborative (VIIM) (ALC). Variant modeling was supported by funding from SciLifeLab's Pandemic Laboratory Preparedness program to B.M. (VC-2022-0028) and by the European Union's Horizon 2020 research and innovation programme under grant agreement no. 101003653 (CoroNAb) to B.M.

Persistence of a Frameshifting Deletion in SARS-CoV-2 ORF7a for the Duration of a Major Outbreak.

Australia experienced widespread COVID-19 outbreaks from infection with the SARS-CoV-2 Delta variant between June 2021 and February 2022. A 17-nucleotide frameshift-inducing deletion in ORF7a rapidly became represented at the consensus level (Delta-ORF7aΔ17del) in most Australian outbreak cases. Studies from early in the COVID-19 pandemic suggest that frameshift-inducing deletions in ORF7a do not persist for long in the population; therefore, Delta-ORF7aΔ17del genomes should have disappeared early in the Australian outbreak. In this study, we conducted a retrospective analysis of global Delta genomes to characterise the dynamics of Delta-ORF7aΔ17del over time, determined the frequency of all ORF7a deletions worldwide, and compared global trends with those of the Australian Delta outbreak. We downloaded all GISAID clade GK Delta genomes and scanned them for deletions in ORF7a. For each deletion we identified, we characterised its frequency, the number of countries it was found in, and how long it persisted. Of the 4,018,216 Delta genomes identified globally, 134,751 (~3.35%) possessed an ORF7a deletion, and ORF7aΔ17del was the most common. ORF7aΔ17del was the sole deletion in 28,014 genomes, of which 27,912 (~99.6%) originated from the Australian outbreak. During the outbreak, ~87% of genomes were Delta-ORF7aΔ17del, and genomes with this deletion were sampled until the outbreak's end. These data demonstrate that, contrary to suggestions early in the COVID-19 pandemic, genomes with frameshifting deletions in ORF7a can persist over long time periods. We suggest that the proliferation of Delta-ORF7aΔ17del genomes was likely a chance founder effect. Nonetheless, the frequency of ORF7a deletions in SARS-CoV-2 genomes worldwide suggests they might have some benefit for virus transmission.

Differential transcriptomic responses to heat stress in surface and subterranean diving beetles.

Subterranean habitats are generally very stable environments, and as such evolutionary transitions of organisms from surface to subterranean lifestyles may cause considerable shifts in physiology, particularly with respect to thermal tolerance. In this study we compared responses to heat shock at the molecular level in a geographically widespread, surface-dwelling water beetle to a congeneric subterranean species restricted to a single aquifer (Dytiscidae: Hydroporinae). The obligate subterranean beetle Paroster macrosturtensis is known to have a lower thermal tolerance compared to surface lineages (CTmax 38 °C cf. 42-46 °C), but the genetic basis of this physiological difference has not been characterized. We experimentally manipulated the thermal environment of 24 individuals to demonstrate that both species can mount a heat shock response at high temperatures (35 °C), as determined by comparative transcriptomics. However, genes involved in these responses differ between species and a far greater number were differentially expressed in the surface taxon, suggesting it can mount a more robust heat shock response; these data may underpin its higher thermal tolerance compared to subterranean relatives. In contrast, the subterranean species examined not only differentially expressed fewer genes in response to increasing temperatures, but also in the presence of the experimental setup employed here alone. Our results suggest P. macrosturtensis may be comparatively poorly equipped to respond to both thermally induced stress and environmental disturbances more broadly. The molecular findings presented here have conservation implications for P. macrosturtensis and contribute to a growing narrative concerning weakened thermal tolerances in obligate subterranean organisms at the molecular level.

SARS-CoV-2 Omicron BA.5: Evolving tropism and evasion of potent humoral responses and resistance to clinical immunotherapeutics relative to viral variants of concern.

BACKGROUND: Genetically distinct viral variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have been recorded since January 2020. The introduction of global vaccine programs has contributed to lower COVID-19 hospitalisation and mortality rates, particularly in developed countries. In late 2021, Omicron BA.1 emerged, with substantially altered genetic differences and clinical effects from other variants of concern. Shortly after dominating global spread in early 2022, BA.1 was supplanted by the genetically distinct Omicron lineage BA.2. A sub-lineage of BA.2, designated BA.5, presently has an outgrowth advantage over BA.2 and other BA.2 sub-lineages. Here we study the neutralisation of Omicron BA.1, BA.2 and BA.5 and pre-Omicron variants using a range of vaccine and convalescent sera and therapeutic monoclonal antibodies using a live virus neutralisation assay. Using primary nasopharyngeal swabs, we also tested the relative fitness of BA.5 compared to pre-Omicron and Omicron viral lineages in their ability to use the ACE2-TMPRSS2 pathway. METHODS: Using low passage clinical isolates of Clade A.2.2, Beta, Delta, BA.1, BA.2 and BA.5, we determined humoral neutralisation in vitro in vaccinated and convalescent cohorts, using concentrated human IgG pooled from thousands of plasma donors, and licensed monoclonal antibody therapies. We then determined infectivity to particle ratios in primary nasopharyngeal samples and expanded low passage isolates in a genetically engineered ACE2/TMPRSS2 cell line in the presence and absence of the TMPRSS2 inhibitor Nafamostat. FINDINGS: Peak responses to 3 doses of BNT162b2 vaccine were associated with a 9-fold reduction in neutralisation for Omicron lineages BA.1, BA.2 and BA.5. Concentrated pooled human IgG from convalescent and vaccinated donors and BNT162b2 vaccination with BA.1 breakthrough infections were associated with greater breadth of neutralisation, although the potency was still reduced 7-fold across all Omicron lineages. Testing of clinical grade antibodies revealed a 14.3-fold reduction using Evusheld and 16.8-fold reduction using Sotrovimab for the BA.5. Whilst the infectivity of BA.1 and BA.2 was attenuated in ACE2/TMPRSS2 entry, BA.5 was observed to be equivalent to that of an early 2020 circulating clade and had greater sensitivity to the TMPRSS2 inhibitor Nafamostat. INTERPRETATION: Observations support all Omicron variants to significantly escape neutralising antibodies across a range of vaccination and/or convalescent responses. Potency of therapeutic monoclonal antibodies is also reduced and differs across Omicron lineages. The key difference of BA.5 from other Omicron sub-variants is the reversion in tropism back to using the well-known ACE2-TMPRSS2 pathway, utilised efficiently by pre-Omicron lineages. Monitoring if these changes influence transmission and/or disease severity will be key for ongoing tracking and management of Omicron waves globally. FUNDING: This work was primarily supported by Australian Medical Foundation research grants MRF2005760 (ST, GM & WDR), MRF2001684 (ADK and ST) and Medical Research Future Fund Antiviral Development Call grant (WDR), Medical Research Future Fund COVID-19 grant (MRFF2001684, ADK & SGT) and the New South Wales Health COVID-19 Research Grants Round 2 (SGT).

Detection of Ganciclovir-Resistant Cytomegalovirus in a Prospective Cohort of Kidney Transplant Recipients Receiving Subtherapeutic Valganciclovir Prophylaxis.

Cytomegalovirus infection during antiviral prophylaxis occurs in transplant recipients despite individualized regimens based on renal function. Fifty kidney transplant recipients were assessed between 2016 and 2019 for valganciclovir dosing, ganciclovir exposure, cytomegalovirus infection, and genotypic resistance markers during the first year posttransplant. Ganciclovir plasma concentrations were measured using mass spectrometry. Population pharmacokinetics was used to determine individual ganciclovir exposure and to evaluate the ability of manufacturer dosing guidelines to meet therapeutic target daily area under the curve (AUC24) of 40 to 50 µg·h/mL. Full-length UL54 and UL97 were assessed using high-throughput sequencing in cytomegalovirus DNA-positive patient specimens. Valganciclovir doses administered to recipients with creatinine clearance of <40 mL/min were higher than specified by guidelines, and they were lower for recipients with creatinine clearance of ≥40 mL/min. The mean ganciclovir AUC24 was 33 ± 13 µg·h/mL, and 82% of subjects did not attain the therapeutic target. Pharmacokinetic simulations showed that the guidelines similarly could not attain the therapeutic target in 79% of individuals. Cytomegalovirus breakthrough occurred in 6% (3/50) of recipients, while 12% (6/50) developed late-onset infection. The mean AUC24s of recipients with (n = 3) and without (n = 47) infection were not significantly different (P = 0.528). However, one recipient with an AUC24 of 20 µg·h/mL acquired two UL97 ganciclovir resistance mutations. Current prophylaxis guidelines resulted in subtherapeutic ganciclovir exposure in several study recipients, including the emergence of resistance genotypes. IMPORTANCE This study examined the pharmacokinetics and viral genomic data from a prospective cohort of kidney transplant recipients undergoing valganciclovir prophylaxis for cytomegalovirus (CMV) prevention. We showed for the first time using high-throughput sequencing the detection of ganciclovir resistance mutations in breakthrough CMV infection during subtherapeutic plasma ganciclovir as indicated by the pharmacokinetic parameter daily area under the curve (AUC24). In addition, we found that current valganciclovir dosing guidelines for CMV prophylaxis are predicted to attain therapeutic targets in only 21% of recipients, which is consistent with previous pharmacokinetic studies. The novel findings of resistance mutations during subtherapeutic ganciclovir exposure presented here can inform future studies investigating the dynamics of drug selection pressure and the emergence of resistance mutations in vivo.

Platform for isolation and characterization of SARS-CoV-2 variants enables rapid characterization of Omicron in Australia.

Genetically distinct variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have emerged since the start of the COVID-19 pandemic. Over this period, we developed a rapid platform (R-20) for viral isolation and characterization using primary remnant diagnostic swabs. This, combined with quarantine testing and genomics surveillance, enabled the rapid isolation and characterization of all major SARS-CoV-2 variants circulating in Australia in 2021. Our platform facilitated viral variant isolation, rapid resolution of variant fitness using nasopharyngeal swabs and ranking of evasion of neutralizing antibodies. In late 2021, variant of concern Omicron (B1.1.529) emerged. Using our platform, we detected and characterized SARS-CoV-2 VOC Omicron. We show that Omicron effectively evades neutralization antibodies and has a different entry route that is TMPRSS2-independent. Our low-cost platform is available to all and can detect all variants of SARS-CoV-2 studied so far, with the main limitation being that our platform still requires appropriate biocontainment.

Different Genes are Recruited During Convergent Evolution of Pregnancy and the Placenta.

The repeated evolution of the same traits in distantly related groups (convergent evolution) raises a key question in evolutionary biology: do the same genes underpin convergent phenotypes? Here, we explore one such trait, viviparity (live birth), which, qualitative studies suggest, may indeed have evolved via genetic convergence. There are >150 independent origins of live birth in vertebrates, providing a uniquely powerful system to test the mechanisms underpinning convergence in morphology, physiology, and/or gene recruitment during pregnancy. We compared transcriptomic data from eight vertebrates (lizards, mammals, sharks) that gestate embryos within the uterus. Since many previous studies detected qualitative similarities in gene use during independent origins of pregnancy, we expected to find significant overlap in gene use in viviparous taxa. However, we found no more overlap in uterine gene expression associated with viviparity than we would expect by chance alone. Each viviparous lineage exhibits the same core set of uterine physiological functions. Yet, contrary to prevailing assumptions about this trait, we find that none of the same genes are differentially expressed in all viviparous lineages, or even in all viviparous amniote lineages. Therefore, across distantly related vertebrates, different genes have been recruited to support the morphological and physiological changes required for successful pregnancy. We conclude that redundancies in gene function have enabled the repeated evolution of viviparity through recruitment of different genes from genomic "toolboxes", which are uniquely constrained by the ancestries of each lineage.

Assessment of Inter-Laboratory Differences in SARS-CoV-2 Consensus Genome Assemblies between Public Health Laboratories in Australia.

Whole-genome sequencing of viral isolates is critical for informing transmission patterns and for the ongoing evolution of pathogens, especially during a pandemic. However, when genomes have low variability in the early stages of a pandemic, the impact of technical and/or sequencing errors increases. We quantitatively assessed inter-laboratory differences in consensus genome assemblies of 72 matched SARS-CoV-2-positive specimens sequenced at different laboratories in Sydney, Australia. Raw sequence data were assembled using two different bioinformatics pipelines in parallel, and resulting consensus genomes were compared to detect laboratory-specific differences. Matched genome sequences were predominantly concordant, with a median pairwise identity of 99.997%. Identified differences were predominantly driven by ambiguous site content. Ignoring these produced differences in only 2.3% (5/216) of pairwise comparisons, each differing by a single nucleotide. Matched samples were assigned the same Pango lineage in 98.2% (212/216) of pairwise comparisons, and were mostly assigned to the same phylogenetic clade. However, epidemiological inference based only on single nucleotide variant distances may lead to significant differences in the number of defined clusters if variant allele frequency thresholds for consensus genome generation differ between laboratories. These results underscore the need for a unified, best-practices approach to bioinformatics between laboratories working on a common outbreak problem.

Simple and reliable protocol for identifying talented junior players in team sports using small-sided games.

We designed and tested a protocol for measuring the performance of individuals in small-sided soccer games. We tested our protocol on three different groups of youth players from elite Brazilian football academies. Players in each group played a series of 3v3 games, in which individuals were randomly assigned into new teams and against new opponents for each game. We calculated each individual's average individual goals scored, goals scored by teammates, goals conceded, and net team goals per game. Our protocol was consistent across days and repeatable across groups, with intraclass correlation coefficients (ICCs) of 0.57-0.69 for average net goals per game across testing days. Players could achieve high success by scoring goals or ensuring their team concede few goals. We also calculated the first and second dimension of a principal component analysis based on each player's number of goals scored, goals scored by teammates, and number of goals conceded per game. Players that were overall high performers had higher PC1 scores, while PC2 scores represented the type of contribution made by a player to overall performance. Positive PC2 values were indicative of high number of individual goals while negative values were associated with more goals from teammates and fewer conceded goals. Our design allows coaches and scouts to easily collect a robust metric of individual performance using randomly designed, small-sided games. We also provide simulations that allow one to apply our methodology for individual talent identification to other small-sided games in any team sport.

Phylogenomic Insights into Deep Phylogeny of Angiosperms Based on Broad Nuclear Gene Sampling.

Angiosperms (flowering plants) are the most diverse and species-rich group of plants. The vast majority (∼99.95%) of angiosperms form a clade called Mesangiospermae, which is subdivided into five major groups: eudicots, monocots, magnoliids, Chloranthales, and Ceratophyllales. The relationships among these Mesangiospermae groups have been the subject of long debate. In this study, we assembled a phylogenomic dataset of 1594 genes from 151 angiosperm taxa, including representatives of all five lineages, to investigate the phylogeny of major angiosperm lineages under both coalescent- and concatenation-based methods. We dissected the phylogenetic signal and found that more than half of the genes lack phylogenetic information for the backbone of angiosperm phylogeny. We further removed the genes with weak phylogenetic signal and showed that eudicots, Ceratophyllales, and Chloranthales form a clade, with magnoliids and monocots being the next successive sister lineages. Similar frequencies of gene tree conflict are suggestive of incomplete lineage sorting along the backbone of the angiosperm phylogeny. Our analyses suggest that a fully bifurcating species tree may not be the best way to represent the early radiation of angiosperms. Meanwhile, we inferred that the crown-group angiosperms originated approximately between 255.1 and 222.2 million years ago, and Mesangiospermae diversified into the five extant groups in a short time span (∼27 million years) at the Early to Late Jurassic.

Paternally-biased gene expression follows kin-selected predictions in female honey bee embryos.

The Kinship Theory of Genomic Imprinting (KTGI) posits that, in species where females mate with multiple males, there is selection for a male to enhance the reproductive success of his offspring at the expense of other males and his mating partner. Reciprocal crosses between honey bee subspecies show parent-of-origin effects for reproductive traits, suggesting that males modify the expression of genes related to female function in their female offspring. This effect is likely to be greater in the Cape honey bee (Apis mellifera capensis), because a male's daughters have the unique ability to produce female offspring that can develop into reproductive workers or the next queen without mating. We generated reciprocal crosses between Capensis and another subspecies and used RNA-seq to identify transcripts that are over- or underexpressed in the embryos, depending on the parental origin of the gene. As predicted, 21 genes showed expression bias towards the Capensis father's allele in colonies with a Capensis father, with no such bias in the reciprocal cross. A further six genes showed a consistent bias towards expression of the father's allele across all eight colonies examined, regardless of the direction of the cross. Consistent with predictions of the KTGI, six of the 21 genes are associated with female reproduction. No gene consistently showed overexpression of the maternal allele.

Emergence of an evolutionary innovation: Gene expression differences associated with the transition between oviparity and viviparity.

Our understanding of the evolution of complex biological traits is greatly advanced by examining taxa with intermediate phenotypes. The transition from oviparity (egg-laying) to viviparity (live-bearing) has occurred independently in many animal lineages, but there are few phenotypic intermediates. The lizard Saiphos equalis exhibits bimodal reproduction, with some viviparous populations, and other oviparous populations with long egg-retention, a rare trait where most of embryonic development occurs inside the mother prior to late ovipositioning. We posit that oviparous S. equalis represent an intermediate form between "true" oviparity and viviparity. We used transcriptomics to compare uterine gene expression in these two phenotypes, and provide a molecular model for the genetic control and evolution of reproductive mode. Many genes are differentially expressed throughout the reproductive cycle of both phenotypes, which have clearly different gene expression profiles overall. The differentially expressed genes within oviparous and viviparous individuals have broadly similar biological functions putatively important for sustaining embryos, including uterine remodelling, respiratory gas and water exchange, and immune regulation. These functional similarities indicate either that long egg-retention is an exaptation for viviparity, or might reflect parallel evolution of similar gravidity-related changes in gene expression in long egg-retention oviparity. In contrast, gene expression changes across the reproductive cycle of long egg-retaining oviparous S. equalis are dramatically different from those of "true" oviparous skinks (such as Lampropholis guichenoti), supporting our assertion that oviparous S. equalis exhibit an intermediate phenotype between "true" oviparity and viviparity.

Accounting for Uncertainty in the Evolutionary Timescale of Green Plants Through Clock-Partitioning and Fossil Calibration Strategies.

Establishing an accurate evolutionary timescale for green plants (Viridiplantae) is essential to understanding their interaction and coevolution with the Earth's climate and the many organisms that rely on green plants. Despite being the focus of numerous studies, the timing of the origin of green plants and the divergence of major clades within this group remain highly controversial. Here, we infer the evolutionary timescale of green plants by analyzing 81 protein-coding genes from 99 chloroplast genomes, using a core set of 21 fossil calibrations. We test the sensitivity of our divergence-time estimates to various components of Bayesian molecular dating, including the tree topology, clock models, clock-partitioning schemes, rate priors, and fossil calibrations. We find that the choice of clock model affects date estimation and that the independent-rates model provides a better fit to the data than the autocorrelated-rates model. Varying the rate prior and tree topology had little impact on age estimates, with far greater differences observed among calibration choices and clock-partitioning schemes. Our analyses yield date estimates ranging from the Paleoproterozoic to Mesoproterozoic for crown-group green plants, and from the Ediacaran to Middle Ordovician for crown-group land plants. We present divergence-time estimates of the major groups of green plants that take into account various sources of uncertainty. Our proposed timeline lays the foundation for further investigations into how green plants shaped the global climate and ecosystems, and how embryophytes became dominant in terrestrial environments.